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Abstract: We study out-of-equilibrium quasi-particle distributions of the Higgs 
and W fields during tlie zero-temperature tachyonic electroweak transition that has 
been assumed in recent scenarios of baryogenesis. Approximating the process by 
a fast quench, we perform classical real-time lattice simulations in the SU(2)-Higgs 
model. The emerging quasi-particle numbers and energies are then used to determine 
the effective temperatures, chemical potentials and masses of the particles shortly 
after the transition. 
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1. Introduction 

An electroweak transition, in which the particles of the Standard Model acquired 
their masses, is assumed to have taken place in the early universe. According to the 
standard lore, it was a finite-temperature transition. However, in recent scenarios 
of electroweak baryogenesis |1|, ^, this transition is assumed to have taken place at 
essentially zero temperature shortly after low-scale inflation ||^ by the effective mass- 
squared parameter of the Higgs field going negative ('tachyonic') ^ ^ 0]. The 
resulting spinodal instability has been shown to provide an effective mechanism of 
preheating i, |, 0,|ll|,|ll. 

These are topics of non-equilibrium field theory. Quantum fields that are way 
out of equilibrium need to be treated non-perturbatively, which is well-known to be 
a difficult task. When applicable, classical approximations can be very useful, since 
they can be treated by numerical simulation. The tachyonic electroweak transition 
is an excellent example, since the classical approximation is well justified 
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Whereas results so obtained are within the language of classical fields, it is desir- 
able to connect with the terminology of kinetic theory, when appropriate. Particle- 
distribution functions have intuitive appeal; they may be describable by time-depend- 
ent effective temperatures and chemical potentials. In this paper we address the 
problem of extracting particle-distribution functions from field-correlation functions 
obtained in classical approximations. 

As is well known, the identification of quasi-particle distributions in field theory 
is not unique. We use a method that was introduced in ||13| for fermions and |T^ for 
bosons. This method has also been found useful in other out-of-equilibrium studies 
[ pTS] , p!6| , \r7\ , [T8|, |19|, ^ and in QCD studies using the classical approximation 
applied to the initial stage of heavy-ion collisions [P2l Bl, 123^ . 



An important topic in electroweak baryogenesis is the time needed for the sys- 
tem to reach approximate thermalization after the transition, and the value of the 
corresponding effective temperature. The thermalization should be fast enough and 
the temperature low enough to prevent a possible washout of the generated baryon 
asymmetry by sphaleron transtions. One of the results in this paper is an estimate 
of this effective-thermalization time and temperature. However, we shall also find 
that a temperature is not sufficient to characterize the particle distribution and that 
a substantial chemical potential is needed as well. A preliminary application to the 
present case is in . 



In section 2 we introduce the equations of motion and recall the initial conditions 
for the tachyonic electroweak quenching transition [^. Section 3 deals with the 
definition of distribution functions for the Higgs and W particles, with more details 
in the Appendix. In section 4 we present results of numerical simulations: particle 
numbers and effective energies and the determination of effective temperatures and 
chemical potentials. A discussion of the results is in section 5. 



2. Tachyonic electroweak transition 



In a tachyonic electroweak transition the effective mass-squared parameter of the 
Higgs field in the effective potential is assumed to change sign from positive to 
negative. In hybrid inflation models 0, 0, this is caused by the coupling of the 
Higgs fleld to the inflaton fleld. In a flrst exploration we assume the transition to be 
dominated by the dynamics of the SU(2) gauge- Higgs sector of the Standard Model, 
and model the transition by a quench. At the electroweak energy scale the expansion 
rate of the universe (0(10^^) eV) is negligible compared to the dynamical time scales 
of the flelds and the process can be studied in Minkowski space-time. 
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2.1 Equations of motion 

The SU(2)-Higgs model is given by the action 

1 



Sr i = - d X 



■^9 



(2.1) 



with F^^ = d^A^-duAf,-i[A^,A^], = (d^-iA^)^. Furthermore, A^ = A;^rV2, 
the r°, a = 1,2,3 are the Pauli matrices and is the Higgs doublet (our metric is 
diag(— 1, 1, 1, 1)). The zero-temperature Higgs and W masses are given by = 
2Af^ = 2/i^, = g'^v'^/A, with v = fi/VX the vacuum expectation value of the 
Higgs field. The equations of motion are (m, n = 1, 2, 3): 

9o'0 = DnDn^ + _ 2A(0V)0, (2.2) 
doE!^ = D^Ft^ + !|![(D„0)tr> - 4>^T^D,,4>1 (2.3) 

with E'^ = F°Q, Z)^* the covariant derivative in the adjoint representation and we have 
chosen the temporal gauge, Aq = 0. The equations of motion for Aq constitute the 
three Gauss-constraint equations, which are to be satisfied by the initial conditions. 
They are conserved by the equations of motion, and read 

D'^tK = ^(^o^V'^^ - <pW''do<P). (2.4) 

For the numerical simulations the action is discretized on a space-time lattice, 
which leads to discretized equations of motion, see e.g. |2^. More details are given 
in 
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2.2 Initial conditions 

Before the electroweak transition the system is assumed to be in the symmetric phase 
{{(j)) = 0) corresponding to an effective action in which the term —^'^(p'^cf) in (|2.1j ) is 
replaced by filfi(j)'^(f), with positive /^^g. The transition is then caused by /Xgg going 
through zero and ending up at today's value — /U^. We model this process by a quench, 
in which /i^g has magnitude /i^ and fiips its sign instantaneously: /Xgg = /x^ ^ ~/^^. 
This approximation gives maximal out-of-equilibrium conditions, which we use for 
testing the baryogenesis scenario |]^, ^ . The more gradual transition expected from 
the coupling to the infiaton field gives qualitatively similar results ||^, |^ . 

The state just after the quench is unstable (the spinodal instability), and in 
the limit X —>■ 0, g'^ it is possible to solve exactly the time evolution of the 
momentum modes of the Higgs field 

cj){k,t) = [ d^xe-'^-''(j){x,t). (2.5) 



/L3 

Here is the volume of our system with periodic boundary conditions 
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The low momentum modes (|k| < /i) initially grow roughly exponentially until 
the interaction terms in the equations of motion become important. Before that 
happens, the exponential growth leads to large occupation numbers and to effectively 
sharp values of both canonical variables and vr = c?o0, which justifies the subsequent 
use of the classical approximation |]28|, ^, Q . The full nonlinear back-reaction is thus 
taken into account without further approximation. This classical approximation 
should be reasonable for observables that are dominated by the low momentum 
modes, until classical equipartition sets in, which, as we shall see, does not occur on 
the time scale of our simulation. 

In it was shown that a consistent way to initialize the Higgs field (in the 
initially free-field approximation), is to generate classical realizations of an ensemble 
that reproduces the quantum vacuum correlators in the symmetric phase before the 
quench, i.e. 



(0(k)0(k)t) = (7r(k)7r(k)t) = y^E+E, (2.6) 

However, we only initialize the unstable (low momentum k < fi) modes {k = |k|). 
This initial condition scheme is the "Just the half" case of 0. 

There is an issue concerning how to choose these initial realizations such that 
they obey the global Gauss law. This technical point is treated in Appendix A of 
our paper [27| and, for the case of 1+1 dimensions, in [^. Given a realization of the 
Higgs field with zero total isospin charges, we can solve the local Gauss-law equations 
(|]D to find the . We set A'} = initially. 

Because only the modes with k < fi are initialized, we do not expect problems 
with lattice artifacts until the system equilibrates classically. We will see that this 
happens on timescales much longer than the ones reached in our simulations (see 



also 25, 



3. Distribution functions 

Under appropriate circumstances kinetic equations can be derived in field theory, 
in which particle numbers constitute a reduced set of dynamical variables, see, e.g. 
[pO[] for the case of QCD. Here we are not concerned with this role of distribution 
functions, since the dynamics is treated numerically, but consider them as observables 
for studying the preheating process after the quench. 

For interacting fields out of equilibrium the definition of local particle numbers 
is not unique. At finite (and typically short) times it is not possible to 'go on shell', 
and there may be damping and finite-width effects. However, the system may display 
effective particle-like behavior in the two-point correlation functions. 

In non-abelian gauge theories there is also the question of gauge invariance when 
the chosen two-point functions are not gauge invariant, as is usually the case. One 
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Figure 1: Higgs field expectation values {cj)^(f))/v^ versus time. Inset: the long time 
behavior. 



may render them gauge invariant by supplying parallel transporters along suitable 
paths, as advocated in |Q . This introduces path-dependence. More generally, there 



may be field dependence: different fields with the same quantum numbers have 
different correlators and the resulting particle numbers may or may not depend on 
the choice being made. 

Part of the ambiguities reside in the identification of effective-particle energies, 
the dispersion relations. We use a method |T^, [1^ in which the effective-particle 
numbers and energies are determined selfconsistently. Ref. |T^ also contains a study 



of the effect of using parallel transporters and makes a comparison with the Wigner- 
function approach. See also |^ for further details on the fermionic case and the 



relation with conserved charges. 

3.1 Higgs and W particle numbers and effective energies 

We are interested in the Higgs- and W-quasi-particle distributions, which can be 
obtained from the 0- and A-correlation functions in a suitable gauge. The natural 
choice of gauge for the Higgs fields is the unitary gauge, in which has only one 
non-zero real component. For the gauge fields, we will study the particle distribution 
in both the unitary gauge and in the Coulomb gauge diAi = 0. 
Writing the Higgs field in the form 



y/2 VV54(x,t) - iv53(x,t)/ 



x,t)=^ rT r ^ (3.1) 
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where 99^, a = 1, 



4 are real, the unitary gauge is defined by 

<^^ = 0, a = 1,2, 3, /i = (^4>0. (3.2) 
The normahzation of the fields is chosen such that, in 



Accordingly, h = ^J2^^). 
the small amplitude approximation around a ground state, the fields enter in the 
kinetic part of the (in general effective) action with the canonical normalization. In 
the present case we have S = J d^x {l/2)dohdoh + ■ ■ ■ . After extraction of the gauge 
coupling g, A'^ ^ Q^ni ^^e field A'^ is also properly normalized. This normalization 
criterion can also be applied to composite fields, e.g. the rho-meson field in QCD, 
provided their effective-action approximation is known. 

For simplicity, let us first neglect the coarse graining that is implicit to the 
notion of distribution functions, and come back to this later. The particle numbers 
and energies are defined [|l^, 0] by analogy with the free-field expectation values 
of the equal-time correlators (00), (tttt), (^^4), {EE) (see also the Appendix). The 



(3.3) 
(3.4) 



Higgs particle numbers and effective energies are defined as 

(t) = v/(/i(k, t)/i(-k, t)),(7r^(k, t)7r;,(-k, t)). 



f(t) 



(7r?,(k,t)7r?,(-k,t)), 



(/i(k,t)/i(-k,t)), ' 

where (■ ■ ■)c is the connected two-point function given by {AB)c = {AB) — {A){B). 
We have replaced the 'quantum' particle- number combination Uk + 1/2 with the 
'classical' n^, and used spherical symmetry to write these as a function of A; = |k|. 

In a general gauge, we can decompose the gauge field two-point function in a 
transverse and a longitudinal part, as 



{At{k,t)A'^{-k,t)) 
{E:ik,t)E'i-k,t)) 



S, 



ab 



6. 



ab 



D^{k,t) + ^Dt{k,t) 
D^{k,t) + ^Df{k,t) 



(3.5) 
(3.6) 



The Sab refiect the fact that the initial conditions are isospin symmetric. In the 
Coulomb gauge, the gauge potential is purely transverse, D^{k, t) = (but D^{k, t) 7^ 
0), and the and cOk can be defined analogously to the Higgs case: 



ntit)^jD^ik,t)D^ik,t) 



lD^{k,t) 
D^{k,t) 



(3.7) 



In the unitary gauge, the free-field correlators are those of a massive Yang-Mills field 
(the derivation is given in the Appendix): 



{A{k,t)Aj{-k,t)) 
{E,{k,t)E,i-k,t)) 



6ij + 


kikj \ Uk 


2 


m^J ujk 


Sij - 


kikj \ 


P + m^y 



UkUJk ■ 



(3.8) 
(3.9) 
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We decompose this into longitudinal and transverse modes, and arrive at the expres- 
sions we will use to determine longitudinal and transverse occupation numbers and 
mode energies in the unitary gauge, 



D^Dm, ^^ = \l^y (3.10) 



nt = \lDt{k)Dl{k) , uj\: = m^l^l^^^. (3.11) 

The transverse case is analogous to the Higgs case. In the longitudinal case we 
assumed the form = m^g + k"^ and replaced m^g in ( pl8| , |3l9| ) . Then ( p.llj ) 

follows straightforwardly. Note the inverse dependence on / in uj^ . In practice 
we analyze the data by first replacing m^^ — ^ mw in ( |3.11| ), and then correct for it. 



3.2 Coarse graining 

Consider the problem of defining a position-dependent distribution function n(x, k, t) 
for a system out of spatial equilibrium. A natural approach is to consider a region 
-R(x) of size B around the position x, e.g. a cube of volume B^, and to focus on 
this region. This means restricting the Fourier integrals etc. in the formulas in the 
previous section to the region i?(x). The size of the region then determines the 
precision in momenta and positions of the particles. Similarly, one expects to have 
to do some coarse graining in time in order to control the fluctuations in the energies 
of the particles. Such time averages were taken in 0, |l6l , but since we try to follow 
the out-of-equilibrium process in time as closely as possible, we shall not do so here. 

For simplicity, consider a scalar field in 1+1 dimensions and let the localized 
region R{x) be the interval {x — B/2,x + B/2). The correlators in momentum space 
associated with R{x) are then given by (suppressing the common time label) 

r-x+B/2 ^—iky+ikz 

C^^{x,k) = / dydz {ip{y)(p{z))c, 

Jx-B/2 -D 
nx+B/2 ^-iky+ikz 

C^^{x,k) = / dydz (7r(y)7r(2;))c, (3.12) 

Jx-B/2 -D 

from which we obtain the distribution functions in the quantum theory as 
n{x,k) + 1/2 
uj{x, k) 

In the classical approximation the '1/2' is left out. 

For a homogeneous system we can improve statistics by taking a spatial average, 

C^^{k) = Y I dxC^^{x,k) = ^w{p- k) {if{p)if{-p)), (3.14) 

p 

4sin^[(p-A;)i?/2] 



k) = ^ ^-(^' = (^(^' + 1/2) uj{x, k). (3.13) 
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and similar for the tttt correlator {(p{p) is the Fourier transform in the total volume, 
as in (|2.5| )). The weight function w{p — k) is sharply peaked about p = k and 
normalized, 

^w{p-k) = l. (3.16) 



The particle numbers are obtained after the spatial averaging, 

c^^{k) = !^h±lll^ c^^^k) = K + 1/2) ^fc. 



(3.17) 



The coarse graining has the effect of smoothing out the correlators in momentum 
space. This is a welcome feature, since the initial momentum modes of the fields are 
uncorrelated, no matter how large the volume (their variance is given by ( p.6|) ). 

In practise we implement spatial coarse graining by 'binning' our momenta spher- 
ically, for example for the Higgs field: 



(Mk,t)M-k,t)) ^ (Mp,^)M-p,t)), 



(3.18) 



-fc-A<|p|<fe+A 



where A^^ is the number of independent momenta in the momentum bin labelled by 
k. In position space this corresponds to spherical shells of thickness of order tt/A. 

3.3 Early time 

The exponential growth of the particle numbers after the quench was studied in P| 
in the approximation \ = g = Q. For each real mode (pa-, a = 1, ■ ■ ■ , 4, of the Higgs 
field, the particle number in the unstable region {k < /x) is given by^ 



1 



+ 



sinh^(2v//i2 - F t) 



nl/2 



1 



4 



(3.19) 



4 4(^4 _fc4) 

where the last form holds for 2 a/ fi"^ — k'^t ^ 1. The field correlators are given by 

1 



5, 



af3 



l + ^^^sinh^(V;:^t) 



(3.20) 



Using the above information we can estimate the particle number at a time tni where 
the neglected nonlinearities may be expected to stop the exponential growth. We 
identify with the time where (v?aV^o) (unstable modes only) reaches the inflexion 
point \j? l^tX of the Higgs potential, i.e. 



(27r)3 3A" 



(3.21) 



^Here we correct an error by a factor 1/2 in eq. (3.15) of j^. 
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For the parameters (A = 1/9) used in our simulation this gives tni = 3.32321/2 ^. 
Then (|3.19| ) gives for the particle number of each real mode at zero momentum 



n^(ini) = 193. 



(3.22) 



The distribution at that time is plotted in figure H, together with the frequency 
[|]. We have not attempted to compute the corresponding quantities for the radial 



100 




0.2 0.4 0.6 



Figure 2: Left: for a real Higgs mode versus k. Right: ujk and ^/\^j?—W\ versus k. 
The units are = m///\/2 and the time is t = 3.32321 = 4.7 m 



H ■ 



mode h = y/ipa'^a, but expect them to be similar. 



4. Results 

We have performed simulations with g = 2/3, A = 1/9, giving rrijj/rri^ = 2, and vol- 
ume = 21^m]f. All results shown will be quoted in units of the zero-temperature 
Higgs mass rriH- 

We used a 60^ lattice, with lattice spacing a = 0.35 m^^. The maximum mo- 
mentum in each direction is vr/a = 9.0 tuh, but reasonably accurate continuum 
behavior is expected to be limited to the region < 1/a = 2.8m//. As is well- 
known, lattice-artifacts are substantially reduced by using the corrected momentum 
k'i = (2/a) sin(aA;j/2). In the following we always use k^ and drop the prime. 

We generated 42 independent realizations of the initial conditions ( |2.6| ), and 
sampled the subsequent time evolution for tmn = 1,2, . . . , 12, 20, 30, 40, 50, 100. 
The Coulomb-gauge fixing has been performed using an overrelaxation algorithm 
pl| , p2| , stopping when L~^^^^ \diA°;{x)\'^ < 10~^^ in lattice units. This precision 



has been chosen to ensure that also the low-momentum modes are tranverse to a 
high degree of accuracy. 



We have averaged nearby momenta as described in section |3.2| , within 'bins' of 
size 2A = 0.05 a^^ = 0.0175 m//. The zero-mode is always in a bin of its own, and 
will be treated separately in that it will not be included in any fit to the data. 
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Figure 3: Higgs-particle numbers function of k. 
4.1 Particle distributions 

In figure ^ we show the particle distribution for the Higgs fields. We see that the 
occupation numbers of the low-momentum modes increases exponentially up to t ~ 
6 — SmJi^, where they start saturating to an uq of about 100. At this point, the 
high-momentum modes rapidly become populated, and after t 20 m^^ the system 
evolves only very slowly, resulting in an approximately exponential distribution. 

In figure § we show the Coulomb-gauge W-particle numbers for all times. We 
see the same qualitative behavior as for the Higgs field, with the occupation numbers 
increasing exponentially up to t ~ 6mjj^, where they start saturating to no «i 15, 
followed by a rapid growth in the high-momentum modes ('sudden up-sweep of the 
tail') and only a slow evolution for t > 20 m^^. The main difference with the Higgs 
case is that the initial particle numbers are much lower — initially all the energy is 
in the Higgs fields. In figure ^ we show the time evolution of the particle numbers 
uq of the zero modes. 

The exponential growth and subsequent slow evolution can also be seen by plot- 
ting the lowest non-zero momentum modes separately as a function of time, as in 
figure ||. 

Figure |^ shows the W-particle distribution in the unitary gauge, with the longitu- 
dinal and transverse modes plotted separately. Note that the longitudinal zero mode 
is not defined. The same qualitative behavior as before can be seen. In this case the 
initial particle numbers are much larger and close to those of the Higgs field, which 
is natural since the Goldstone modes, which are absorbed into the unitary-gauge 
W- fields, are also populated by the initial conditions ( |2.6| ). 
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Figure 4: Coulomb-gauge W-particle numbers as a function of k. 
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Figure 5: Particle numbers of the zero modes as a function of time; H: Higgs; Wt- 
transverse W in unitary gauge; Wc- transverse W in Coulomb gauge; inset: early time. 



Finally, in figure H, we show the particle distribution for all four 'species' (Coulomb- 
gauge W, transverse and longitudinal unitary-gauge W, and Higgs) at the latest time, 
t = 100 m^^. We see, firstly, that all four have an exponential fall-off at large mo- 
menta. Secondly, for k < 0.7 mj^ the occupation number remains significantly above 
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Figure 6: Coulomb-gauge W-particle numbers for the four lowest momentum modes, as 
a function of time. 
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Figure 7: Unitary-gauge transverse (left) and longitudinal (right) W-particle numbers as 
a function of k. 



1, vindicating our use of the classical approximation. Thirdly, the transverse gauge 
field modes have almost exactly the same distribution in the unitary gauge as in 
the Coulomb gauge. The distribution of the longitudinal modes on the other hand 
deviates somewhat from that of the transverse ones. 

4.2 Dispersion relation and effective mass 

In figure ^ we show the dispersion relation, as a function of k"^, for Coulomb-gauge 
W-particles. For t < 20 mjj^ there is no sensible dispersion relation, as can be seen 
from the data for t = 8; while for t > 30 mjj^ it approaches the form cu^ = m^g + ck"^, 
with c ~ 1. The inset shows the dispersion relation for t = 30, 40, 50 and 100 mjj^ in 
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Figure 8: Particle distributions at the latest time, for all four particle 'species' 
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Figure 9: Dispersion relation for W-particles in Coulomb gauge. 



a more restricted momentum range. The data are very well described by a straight 
line, and for t = 40, 50 and 100 they are indistinguishable. It is striking that the 
data turn out to be compatible with a straight line all the way up to k'^ = 40m^, far 
into the region where one would expect lattice artifacts to dominate and the classical 
approximation to break down. 
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Figure 10: W-particle dispersion relation in the unitary gauge, for transverse (left) and 
longitudinal (right) modes. 



In figure TO we show the dispersion relation in the unitary gauge. In this case, 



a particle-like behavior takes considerably longer to emerge than in the Coulomb 
gauge: for the transverse modes the slope is still smaller than 1 (and increasing) at 
the latest time, t = 100 m];^^, while for the longitudinal modes a curvature remains 
for k"^ > mfj. However, the intercepts (effective mass-squared) are quite compatible. 




Figure 11: Dispersion relation for Higgs particles. Also shown is the best fit to a straight 
line for t = 100 m^^. 



In figure ^ we show the dispersion relation for the Higgs particles, for tmn = 8, 
20, 30, 40, 100. Here again, we find that the dispersion relation is stable for t > 
SOmJj^. Evidently, there is still some remnant of the odd-looking dispersion relation 
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— fc^ I during exponential growth at early times, shown in figure ^. Standard 
+A;^-like behavior emerges between tmu = 12 and 20. Note also that the W- 



eff 



dispersion relation at tmu = 8 in figures ^ and shows similar behavior: the various 
field modes appear to adjust to each other locally in momentum space. This can also 
be seen in the particle numbers. Similar 'local momentum space equilibration' has 



been observed and explained in |TB 



We fit the effective energies to the form uo^ = ck"^ + m^^. The uncertainty in the 
dispersion relation is estimated by varying the end-point /Cmax of the fitting range and 



determining the statistical errors by the bootstrap method. In figure |T2| we show the 
result of this procedure for the slope of the Higgs dispersion relation. We find that in 



1.4 



1.2 



0.8 



0.6 



0.4 



4, 



A-f 





-♦ H, 


t= 20 




♦ - 


-♦ H, 


t= 30 






-♦ H, 


t= 40 






♦ H, 


t= 50 




♦ - 


-♦ H, 


t=100 





1.5 



2.5 



Figure 12: Slope of the dispersion relation for the Higgs field as a function of the maximum 
momentum in the fit. 



order to obtain stable values for c and meff, we need to include points up to /c ~ 2 niH 
in our fits, while beyond this point the fit values do not change significantly. This 
was the case for all the fits we performed, and we have thus chosen fcmax = 2mH 
when quoting the fit parameters and statistical errors in table |I|. For the gauge 
fields, the intercepts (effective masses) in the two gauges are compatible, yielding a 
value rries ~ 0.68 mj^, close to the zero-temperature value 0.71 mjj. However, the 
slope, which in the Coulomb gauge is very close to 1, is considerably lower in the 
unitary gauge, although it appears to approach 1 with increasing time. Quasi-particle 
behavior appears to take longer to emerge in this gauge. For the Higgs field, we find 
an effective mass ~ 0.88 mj^ at t = 100 m^^, which is significantly smaller than the 
zero-temperature Higgs mass, although it appears to be increasing with time. The 
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t 


Wc 


Wt 




H 






c 




c 




c 


30 


0.56(2) 


0.94(2) 


0.59(3) 


0.70(5) 


0.93(2) 


0.91(3) 


40 


0.64(1) 


0.974(6) 


0.67(1) 


0.69(3) 


0.79(3) 


1.02(2) 


50 


0.66(1) 


0.970(7) 


0.67(1) 


0.79(1) 


0.84(2) 


1.01(2) 


100 


0.68(1) 


0.972(7) 


0.69(1) 


0.828(6) 


0.89(2) 


0.98(2) 



Table 1: Effective masses (in mn units) and slopes in the dispersion relation. The errors 
are purely statistical. 



slope is found to be consistent with 1 from t = 40 mj^ onwards. The effective Higgs 
mass appears to be increasing with time, which may be expected from the fact that 
Ugjj = (averaged over an oscillation period) is still slowly increasing (figure 1). 

One would expect a component Ves/v ~ 0.87 in both effective Higgs and W mass 
ratios, but the effective W mass somehow seems to have settled from t = AOmJj^ 
onwards. 

Finally we turn to the longitudinal dispersion relation uj^. In plotting the data 
in figure ^ (right) we used mw instead of m^g in eq. ( |3.11| ), and we now consider 
this point. Let u;^ be the frequency defined by ( |3.11D with m^g — > rriw. Then 

— i^es/''^w)^k "^eff = ("^vk/"^')' whcrc m' = u'q. The data for u'^ with 
k'^/m\ < 1 in figure |l^ (right, tmn = 100) can be fitted well by a straight line 
with slope very close to 1, and effective mass m'^ 0.49 m^. It follows that 

m^g ~ mw to a good approximation. 

The different slopes in the dispersion relations may be interpretable by simple 
parameter changes in an effective quasi-particle lagrangian, but we shall not follow 
up on this here. 

4.3 Approximate thermalization, temperature and chemical potential 
We will model the particle distribution with a Bose-Einstein (BE) distribution. 

It may seem strange to use the BE form instead of the classical = T/{uJk — fJ^), 
but this form would not be able to describe the roughly exponential tail of the Uk 
data. We use the BE form simply as a distribution to compare the data with, in 
order to extract effective temperatures and chemical potentials. The BE form may 
be re-expressed as 

1 \ UJk- fJ- 



Bk^\nil + - =^—^ (4.2) 
V rik/ I 

If we have a Bose-Einstein distribution, is a linear function of uj^, and in a plot 
of B vs uj, the inverse temperature can be read off as the slope. We will therefore 
refer to such plots as 'inverse-temperature plots'. 
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We will use two methods to determine effective temperatures and chemical po- 
tentials. The first (method 1) is to perform a straight-line fit to ([4.2|) , using the data 
for Uk and ook- The second (method 2) is to take Uk from the dispersion relation, 
using the fitted values for c and nies determined in the previous section. We then fit 
Uk directly as a function of Uk to the Bose-Einstein distribution (^.1|). 




a),/m„ 

k H 



Figure 13: Inverse-temperature plot for W particles in Coulomb gauge. 



In figure ^ we show B as a function of to for Coulomb-gauge W particles, for 
t = 20, 30, 40, 50, 100 mjl. For tOk < 2.5 — Sitlh the data are compatible with a Bose- 
Einstein distribution, with a rather large chemical potential. At higher energies, 
as shown by the inset, there is some curvature around Uk = 4. In this region the 
particle numbers are very small and it is clearly way beyond the range of validity for 
the classical approximation. It is nevertheless interesting to see that the qualitative 
behavior of the distribution is unchanged as we go from the infrared to the ultraviolet. 
We also see that the distribution changes only very slightly with time, with the 
effective temperature slowly increasing. The same qualitative picture is found also 
in the unitary gauge. 



Figure |1J shows i? as a function of u for the Higgs fields at the latest time. The 
statistical errors in u are larger here than for the W fields. A fit to a straight line 
through the lowest modes — those where > 1 — deviates from the data at higher 
uj, and, as seen in the inset, there is a curvature in the region Uk = 2.5 — 4. 



Figure [T^ shows the fitted temperature T from the straight-line fits, as a function 
of the end-point cjmax of the fitting range. We have chosen cjmax = 2.1, which is below 
the region where shows curvature in the inserts in figures ITB| and in. 
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Figure 14: Inverse-temperature plot for Higgs particles for t = 100 m^^^. The line is the 
best fit to a straight line for iVk < 1.6 m// 
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Figure 15: Effective temperature, from inverse-temperature fits, as a function of fitting 
range. 



The fit values are given in table |^. Since they were obtained with method 1, 
we have added the subscript 1 to Ti and fii. The effective temperature Ti in the 
unitary gauge turns out to be somewhat lower than in the Coulomb gauge: 0.31 
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t 


Wc 


Wt 




H 




Ti 




Ti 




Ti 




30 


0.320(8) 


0.717(12) 


0.24(4) 


0.73(4) 


0.23(2) 


1.07(3) 


40 


0.319(3) 


0.772(8) 


0.268(8) 


0.742(9) 


0.28(3) 


0.94(3) 


50 


0.323(4) 


0.786(9) 


0.298(3) 


0.758(8) 


0.27(2) 


0.97(3) 


100 


0.345(6) 


0.786(9) 


0.309(4) 


0.781(7) 


0.28(2) 


0.95(3) 



Table 2: Effective temperatures and cliemical potentials, from fits to < cj < 2.1mH- 
The quoted errors are a combination of the statistical errors and systematic errors obtained 
by varying the fit range with Wmax between 1.6 and 2.6m//. 

compared to 0.35 mn at the latest time. This can be put down to the slope of the 
dispersion relation being lower, resulting in a smaller effective energy for the same 
mode. The chemical potentials are consistent between Coulomb and unitary gauge, 
giving a value fi ~ 0.78 m//. This is higher than the effective mass ~ 0.68 m// in 
table |1|, which is of course nonsensical, since it would lead to a pole in rik for very 
small k. The effective temperature for the Higgs is similar to that of the W fields, 
although slightly lower ^ 0.28 m//. The chemical potential for the Higgs comes out 
to be ~ 0.95 m//, again higher than the effective mass ~ 0.88 m// in table [^. 




Figure 16: W-particle number at t = 30, 40, 50, 100 m];^^ as a function of the effective 
energy u), together with the Bose-Einstein distribution ( [4.1| ) with parameters from table |2|. 
The dotted lines indicate the variation allowed by the errors in the fitted T and fi. 



Figure |T6| shows the Coulomb-gauge W-particle number as a function of the 
effective energy uj, together with the best inverse-temperature fit to the five lowest- 



lying modes. These data do not appear to be particularly well described by a Bose- 
Einstein distribution, largely due to what appears to be a somewhat erratic behavior 
of the effective energy Uk- Replacing Uk with values taken from the fits in table |1], we 
obtain a much smoother data plot, as shown in figure O. We fit the 5 lowest-lying 



t 


Wc 


Wt 




H 
















30 


0.50(1) 


0.62(1) 


0.33(3) 


0.64(2) 


0.45(2) 


0.97(2) 


40 


0.420(7) 


0.698(11) 


0.31(1) 


0.713(10) 


0.58(6) 


0.83(3) 


50 


0.427(7) 


0.721(9) 


0.335(7) 


0.719(10) 


0.46(2) 


0.88(2) 


100 


0.424(6) 


0.735(9) 


0.370(5) 


0.733(9) 


0.38(1) 


0.90(2) 



Table 3: Effective temperatures and chemical potentials, from fits to the lowest 5 non-zero 
modes (method 2). 
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Figure 17: Higgs and W-particle numbers at t = 100 as a function of the effective 
energy uj taken from the fitted dispersion relation. Also shown is the Bose-Einstein distri- 
bution (|4.1| ) with parameters from table ^ The zero-modes are outside the boundaries of 
the plot (they are given by (0.68,45(2)) for Wc, (0.69,52(2)) for Wt and (0.89, 12(2)) for 
H). 

modes (for which > 0.5) of these data to a Bose-Einstein distribution (method 
2) and give the results in table 0. The fit describes the displayed data well in the 
important region where the classical approximation is supposed to be valid. 

Again we find that the effective temperature for the W particles in the unitary 
gauge is lower than in Coulomb gauge, due to the lower slope of the dispersion 
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relation. On the other hand, the effective Higgs-temperature is falhng from t = 
AOmJj^ onwards and the unitary-gauge W-temperature is rising. A conversion of 
energy in the Higgs to the W may still be going on at the latest time, which is also 
suggested by the behavior of the zero modes in figure ^. 

The chemical potentials both for the Higgs and W particles are lower than those 
obtained using method 1, reflecting the fact that modes with small uj have a larger 
weight. For the Higgs, fi2 is compatible with the effective mass, but for the W it is 
still larger. Note that the zero-modes have not been included in any of the fits. It 
turns out that they have occupation numbers much lower than would be 'predicted' 
by any Bose-Einstein distribution that would simultaneously fit the other 'classical' 
modes. 

This is another indication (in addition to chemical potentials larger than effective 
masses) that the BE fit will overestimate the true distribution for momenta below 
our lowest finite- volume momentum kmin = 2tt/L = 0.30 mn- The flattening of the 
distribution in figure |^ as A; ^ appears to be very resilient. 

5. Conclusions 

We have obtained Higgs- and W-particle distributions and energies after a quenched 
electroweak transition, when the system is still out of equilibrium, for the case 
my/- The particle distributions could be obtained from early times 
t = l/rriH onwards. On the other hand, the effective energies (frequencies) suffered 
much more from fluctuations and they appeared to get conventional forms only after 
t > 20mjj^. Much more statistics would be required to go determine them with rea- 
sonable accuracy at earlier times (where they retain the memory of the instability). 
The Coulomb-gauge W distribution appears to be smoother than the corresponding 
transverse distribution in the unitary gauge, but the two approach each other and 
by the time 100 they are practically indistinguishable. At this stage the gauge 
dependence has practically disappeared; also the longitudinal W modes have settled 
to nearly the same distribution. We have carried out a few simulations up to times 
500 m];^^, which showed that the distributions change very slowly after time 100 m^^. 

The maximum value of the Higgs-particle number turned out to be smaller than 
the analytical estimate in ( |3.22| ), by a factor of 2. This mismatch is evidently due to 
the fact that with the "Just a half" initial conditions, the interactions (including the 
Coulomb interaction generated by the Gauss constraint) are turned on straight away 
from time zero onwards, and that the driving power of the instability is damped by 
sharing energy with the gauge field. At low momenta the particle numbers are still 
large even at t = 100 mjl and there is no reason to doubt the classical approximation 
for these modes. In contrast, the high- momentum modes are still exponentially 
suppressed at this time and there is no sign yet of classical equipartition (n^ oc 1/k 
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for large k) in our simulation. Hence, lattice artifacts are expected to be small in 
the low-momentum part of the distributions. 

Although we know of no good reason to expect a Bose-Einstein distribution, this 
form fitted the data quite well in the important (low momentum) region (figure [T7|) . 
The same is true at intermediate momenta, where the particle numbers are much 
smaller than one, albeit with slightly different temperature and chemical potential. 
This is fortunate, since it allows for a familiar interpretation of the results. The final 
temperatures of the Higgs and W degrees of freedom turned out to be reasonably 
close at time 100 m^^, which indicates that the system is near kinetic equilibrium. In 
contrast, the rather large chemical potentials show that it is still far from chemical 
equilibrium. For very low momenta (lower than our finite-volume simulation was 
able to deal with) the BE fit has to break down because the chemical potentials are 
slightly larger than the effective masses. 

The temperature in the gauge fields at time 100 is relatively low, 0.44 m//, 
which suggest that sphaleron transitions are highly suppressed. However, such a con- 
clusion cannot be drawn from the temperature alone. The large chemical potentials 
correspond to the fact that there is still quite some power in the low momentum 
modes, and sphaleron transitions might still be non-negligible over a long time span. 
However, in the full Standard Model the Higgs and W particles will have decayed in 
a few hundred mj^^ and we do not see a problem with the baryogenesis scenario. 
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A. Free- field correlators in Higgs and Coulomb gauge 

We derive here the free-field correlators used for guidance in the definition of the 
distribution functions in the unitary gauge and the Coulomb gauge. 

In the unitary gauge, expanding around the ground-state configuration h = v, 
= 0, and keeping only terms up to second order in the fields, leads to the effective 
free-field lagrangian 



^free 



(fx 



-dohdoh — -dnhdnh — -m'jj{h — vY 



_|_ Ip" _ F« _ 177,2 Aa Aa ^ 2 Aa 



(A.l) 



where we rescaled A gA and the nonlinear terms in Ff^^, should be dropped. To 
simplify the notation we suppress the common time label t in the following. 
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The canonical conjugate of the Higgs field is 7r/j(x) = SL/Sdohlx.) = doh{x.). 
Let ujk and Uk be defined in terms of the equal-time correlators, assuming (spatial) 
translation and rotation invariance, by 

(/i(k)/i(-k)), = K + l/2)/a;,, (A.2) 
(7r;,(k)7r,,(-k)), = {rik + 1/2) Uk. (A.3) 

Since TCh = doh, the quantity Uk has the interpretation of a local (in time) frequency. 
The meaning of the quantity Uk is elucidated by introducing the annihilation opera- 
tors 

6(k) = -= K(/i(k) - vS^,o) + ^7r,(k)] , (A.4) 

which satisfy the usual commutation relations with the creation operators &^(k). 
Then, using the canonical commutation relations we get 

(6t(k)6(k)) = nfc + '-{7ih{k)h{-k) - 7r,(-k)/i(k))e 

= rife, (A. 5) 

where the last step follows from rotation invariance. So, Uk is the expectation value 
of the number operator b^{'k)b(k). 

For the gauge field the free- field effective action of is just the sum of three actions 
that are identical in form: one for each value of the isospin index a. For simplicity 
we suppress the index a and also the index W on rriw The canonical momenta 
of the gauge fields are given by 11*^ = SL/SdoAo = 0, n„ = 6L/6dQAn = F^n = 
{doAn — dnAo) = —En. The field Aq is not an independent variable, but it follows 
from the time component of the field equation d^F^^"^ —m'^A'^ — 0, i.e. Aq = (9„n„/m^. 
This can be used to express n„ in terms of do An, 

Smn - n„ = doA„,. (A.6) 

In Fourier space, 

Smn + nn(k) = doAn,{k), n„(k) = ( S^n ' ^O^n(k) (A.7) 

The annihilation operators can be defined by analogy to the scalar case, 

a„(k) = KA„(k) + i9oA„(k)] , (A.8) 

and in this case the canonical commutation relations imply 

K(k), 4 (k')] = i^mn + ^) 4k', (A.9) 
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In terms of spin-polarization vectors e^(k, A), A = 1, 2, 3, satisfying' 

h k 



e;(k, A) ( 5mn - ) e„(k, A') = 5,v, (A.IO) 



^e™(k,A)e:(k,A) = {^5^^ + ^^ , (A.ll) 



A 

the usual annihilation operators for a specific spin state are given by 

a(k, A) = e;(k, A) - «n(k), (A.12) 

a„(k) = J]a(k,A)e„(k,A), (A.13) 

A 

with standard normalization 

[a(k,A),at(k',A')] =5k,k'<5AA'. (A-M) 
In a translation- and rotation-invariant state with 

(a^(k, A)a(k, A)) = n^, independent of A, (A. 15) 
and (am(k)a„(— k) + a,|„(— k)a|j(k)) = 0, we have 

(A^(k)A„(-k)) = (5^r. + {uk + 1/2)-, (A.16) 

(n^(k)n„(-k)) = - -^^^^^ K + 1/2) ^fc, (A.17) 

with ujk = + rn^] recall m = mw in this derivation. 

Conversely, in a more general setting we can follow the same route as for the 
Higgs case, define Uk and cok by ( |A.16| , [A?T7D and then derive ( [A.15| ), which justifies 
the interpretation in terms of particle numbers. Evidently, using rotation invariance 
makes it possible to avoid having to make the assumption that (am(k)a„(— k) + 
a]„(— k)ajj(k)) vanishes. 

Next we consider the Coulomb gauge dnA^ = 0. For this case the free lagrangian 
is given by 

l;dohdoh + Irdoipadoipa - \dnhdnh - - \m]j{h^ - V^) 



^free 



+ \doAldoAl - U^Ald^Al - \r<AlAl 



(A.18) 



^ These are the usual orthonormaUty and completeness relations of a massive vector field with 
mass TO. A realization is given by e„(k, A) = unit vector perpendicular to k, A = 1,2, e„(k, 3) = 
knijJk/km. The vectors can be completed with a time component e°(k. A) = 0, A = 1, 2, e°(k, 3) — 
k/m, such that they satisfy the standard relations e^(k, A)*e^(k, A') — 5x\i, e'^(k, A)e''(k, A) = 
■qf"" + kf'k''/m,'^ and fc^e^(k. A) = 0. 



- 24 - 



where we used dnA'^ = and the parametrization (|3.1| ) of the Higgs field, (f^ = h. 
The canonical momentum of is U'^ = Fq^ = OqA^ — dnAg = —E!^. Again the Aq 
are not independent degrees of freedom but given by the constraint equation 

+ K/K - ^WT^a = 0, (A. 19) 

where tt^ = d^ipa and (9„n^ = — (9„9„ylg. So the Ha mix with the longitudinal com- 
ponents = —dnA^. The ipa degrees of freedom play the role of the longitudinal 
components of the W field, and they can be shown to have a mass my/. We shall not 
give further details as we study those components numerically only in the unitary 
gauge where they are included in the gauge field. 

The canonical commutation relations of transverse components of the gauge field 
have the usual form familiar from Coulomb-gauge electrodynamics: 11'^ = OqA'^, 

[^™(k), nf (-k')] = 6at [Srr^n ' ^) ^kk', (A.20) 

and the equal-time commutator of A'^ with U'^^ equals zero. In an isospin- symmetric, 
translation- and rotation-invariant state we can now define and Uk by 



{Alik)Aii-k)) = 5^, ^5^^ - ^) ^^^±^, (A.21) 

h k 



(nf (k)nf (-k)) = 5,, - (n, + 1/2) (A.22) 

where the frequencies ujk = \/m^ -|- fc^ in the free case. With the usual transverse 
polarization vectors en(k. A), A = 1,2, and annihilation operators 

a(k. A) = e;(k, A)-^ {u^Al + ^^,Al) , (A.23) 

it follows that the have the particle number interpretation 

(at(k,A)a(k,A)) =nfc. (A.24) 
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